Density of states and electron concentration of double heterojunctions subjected to an 

in-plane magnetic field. 



O 

o 
o 

(N 

< 

00 



c/2 



CD. Simserides 

INFM and Scuola Normale Superiore, Piazza dei Cavalieri 7, 1-56126 Pisa, Italia 

We calculate the electronic states of Alj^Gai-jjAs/GaAs/Al^^Gai-i^As double heterojunctions 
subjected to a magnetic field parallel to the quasi two-dimensional electron gas. We study the 
energy dispersion curves, the density of states, the electron concentration and the distribution of 
the electrons in the subbands. 

The parallel magnetic field induces severe changes in the density of states, which are of crucial 
importance for the explanation of the magnetoconductivity in these structures. However, to our 
knowledge, there is no systematic study of the density of states under these circumstances. We 
attempt a contribution in this direction. 

For symmetric heterostructures, the depopulation of the higher subbands, the transition from a 
single to a bilayer electron system and the domination of the bulk Landau levels in the centre the 
wide quantum well, as the magnetic field is continuously increased, are presented in the "energy 
dispersion picture" as well as in the "electron concentration picture" and in the "density of states 
picture" . 
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I. INTRODUCTION 

Althouth the behaviour of the Quasi Two-Dimensional 
Electron Gas (Q2DEG) in the presence of a perpendicu- 
lar magnetic field has been studied extensively, much less 
attention has been devoted to the situation where the 
magnetic field is applied parallel to the Q2DEG. In the 
former case, interesting phenomena e.g. the Shubnikov- 
de Haas effect and the integer ||^ and the fractional 
quantum Hall effects have been observed. In the lat- 
ter case, electrons move under the competing infiuence of 
the Lorenz force and the force due to the quantum well 
confining potential. 

In the presence of an in-plane magnetic field, i?, sin- 
gle heterojunctions [|-||, single @,||, double and 
triple ||l3|] square quantum wells, almost square quan- 
tum wells JT^p^ , asymmetric square quantum wells p6| , 
symmetrical wide single quantum wells |jl^,0 and su- 
perlattices have been considered. 

The experimental studies include single heterojunc- 
tions double square quantum wells [p^|-p^, triple 
square quantum wells [ p3| , wide single quantum wells p^ ] 
and superlatticcs pjj]. The most important experimen- 
tal finding p0| , p^ , ^3|p[ is, according to our opinion, the 
strong conductance "oscillations" due to the variation of 
the density of states (DOS), as B is increased. Conduc- 
tance maxima are identified with depopulations of local 
energy dispersion minima, while conductance minima are 
identified with van Hove singularities at the chemical po- 
tential. This situation has been encountered in symmet- 
rical square double |10 1^ and triple quantum wells 
and in symmetrical wide single quantum wells . While 
in the cases of square double and triple quantum wells 
a simple Tight Binding calculation gave the position of 
the maxima and the minima, a self-consistent calculation 
was indispensable in the case of symmetrical wide single 
quantum wells. 

Theoretical studies of the electronic states are usually 



restricted to simple analytically solvable potential wells, 
to Tight Binding Approximation, or to perturbative ap- 
proximations. Self-consistent studies are up to now a few, 
regarding single heterojunctions thin single quan- 

tum wells jl^ and symmetrical wide single quantum wells 
0- 

In the present work, we use self-consistent calculations 
to study Alj;Gai_2:As/GaAs/Ala;Gai_2:As wide double 
heterojunctions (i.e. a system of two heterojunctions 
with relatively large distance between the two interfaces) 
subjected to an in-plane magnetic field. Below, we sum- 
marize the particular aims of this work. 

Our first aim is to study the density of states when 
the Q2DEG is subjected to an in-plane magnetic field. 
In this case, the DOS is not a step-like function, as it is 
with B=0. We show that its form undergoes important 
changes as B is increased, especially in wide double het- 
erojunctions where usually many subbands are present 
pO| , pT| . The self-consistent study of the electronic states 
and specifically of the DOS is of great importance for 
the explanation of the experimental magnetoconductiv- 
ity in these structures. However, up to now, there is no 
systematic study of the DOS under these circumstances. 
We attempt to give a contribution in this direction. 

Our second aim is to study a bilayer electron system, 
different from the commonly used symmetrical double 
square well. Another potentially bilayer electron system 
is the symmetrical double heterojunctions, when the well 
width is increased a lot, due to the transition from a "per- 
fect" square quantum well to a system of two separated 
heterojunctions |2^. In the former structure a high bar- 
rier separates the two electron layers. In the latter struc- 
trure the barrier is formed from the redistribution of the 
carriers in the well and it is relatively weak. Moreover, 
Smrcka and Jungwirth JItI have shown, by calculating 
the energy dispersion curves in a two-subband situation, 
that symmetrical wide single quantum wells can be po- 
tentially bilayer electron systems when the parallel mag- 
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netic field is increased. Here, we present the depopu- 
lations of the higher subbands and the transition from 
a single layer to a bilayer electron system not only in 
the "energy dispersion picture" , but also in the "electron 
concentration picture". Thus, we calculate the electron 
concentration, n{z) and the distribution of the electrons 
in the subbands, ni{z). We also give the "the density 
of states" picture, which is important for the interpreta- 
tion of the transport experiments. Moreover, we show in 
these three pictures that in the centre of our wide quan- 
tum well, as the magnetic field is further increased, the 
bulk Landau levels dominate. Finally, we give an exam- 
ple of an asymmetric heterostructure. 

The basic theory is presented in section 2 to- 
gether with some analogies to the classical picture. 
In section 3 we present the theoretical results for 
Al2:Gai_2;As/GaAs/Al2;Gai_2;As double heterojunctions 
and we comment on some interesting features observed. 
Our conclusions are summarized in section 4. 



II. BASIC THEORY 



is the electron mass, a is the spin and B is the parallel 
magnetic field. 

U{z) = Ubandoffsetiz) + Uc{z) + Uxc{z) + Ue{z), 

(2.2) 

where Uband offset{z) is the potential energy term due to 
the conduction band minima discontinuity, Uc{z) is the 
Coulombic potential energy, Uxc{z) is the exchange and 
correlation potential energy |Q and Ue{z) is the poten- 
tial energy due to an electric field applied in the z-axis 
e.g. due to a gate. The magnetic field is applied in the 
y-axis i.e. B — (0,5, 0). For the vector potential we 
choose A — {Bz,Q^Q) ||2^. The Hamiltonian becomes: 

Htot = {Px - qBzf/{2m*) + p^/(2m*) -I- pi/ {2m*) + U{z) + 

(2.3) 

We split the spatial and the spin part. ^'(r, (?) = 
il^{r) a{a) and Etot — E^yz ± ^9* IJ-bB. For the spatial 
part the envelope function equation writes: 



When a magnetic field, B, parallel to the y-axis, is 
applied to a three dimensional electron gas, the motion 
in the xz-plane is quantized to Landau levels with en- 
ergy eigenvalues E^z = hoj{i + i), where z is a discrete 
quantum number, h is the reduced Planck constant and 
iv = is the cyclotron frequency, m* is the effective 
mass and q=-e is the electron charge. If we additionally 
apply an electric field, E, in the z-axis, E^z depend not 
only on i, but also on the wavevector in the x-axis, k^. 
Specifically, E^z = nu!{i + ^)- ^i^)"^ -nk^{§). In this 
work we are interested in the configuration with a quan- 
tum well in the z-axis (with or without the electric field in 
the z-axis) and the magnetic field in the y-axis. Again, as 
we discuss below, E^z depend on both i and k^. However, 
generally in this case E^z = Ei(kx) cannot be expressed 
analytically and have to be determined sclf-consistently. 

Of course, without a magnetic field, E^^z — Ei + , 
where now i is the subband index. In all the situations 
described above, the energy eigenvalue in the y-axis is 

E,, — t; — f , where fc„ is the wavevector in the y-axis. 
The spin part of the eigenenergy is E^pin — ±^g* i^bB, 
where g* is the effective Lande factor and /is is the Bohr 
magneton. 

Summarizing, in the present configuration, there is a 
magnetic field in the y-axis, a quantum well in the z-axis 
and possibly an electric field in the z-axis (e.g. an ex- 
ternal field due to a gate). With our choice of axes, the 
Hamiltonian is: 



Htot = [p- qA)yi2mr) + U{z) + g*fiB<f ■ B, (2.1) 

where p is the momentum operator, A is the vector po- 
tential, m* = 0.067TOe is the GaAs effective mass, nie 



[ip, - qBzY/{2m*)+pl/{2m*) +pll{2m*) + U{z)]^ = E^yz^. 

(2.4) 

[px,H] = [py, H] ~ 0. Thus, we look for solutions in the 
form ^j) = -^(^{z)e^^''^ e^^^"^ , where S=hxljy is the area 
of the heterostructure in the xy-plane. The coordinate y 

splits from the coordinates x and z. We have Ey = -2^^ 
while in the xz plane: 



lOz) 
dz^ 



2m 



,eB 



hk 



+ -r[E.. - -^{—Y{z + -^Y U{z)]C{z) = 



eB 



(2.5) 



The non- magnetic part of the potential energy is U(z), 
while the magnetic part of the potential energy is 
^{^Y{z + ^^Y- The center of the magnetic potential 



energy is the point zq = — ^ — Thus, the elec- 

tron is free in the y-axis, but the magnetic field correlates 
the motion in the x-axis and the z-axis. The motion in 
the xz-plane is characterised by a running wave e*'^^^ and 
the bound state Ci,fe^ {z) which depends on both i and kx- 
The energy eigenvalues are: 



1 h^k"^ 1 

Etot = Exz + Ey± -g*fiBB = E,{kx) + tt^ ± xsVb^, 

2 2m 2 



(2.6) 



where, generally Ei{kx) ^ Ei{-kx). 
The density of states is: 
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where: 



2m* 
(2.8) 



We have used the symbohsm E'i^fc^ = Ei{kx)- Intergrat- 
ing over ky, Eq. 8 writes: 



, , LyV2m* 1 



(2.9) 



where O is the step function. We must note here that 
the DOS is not a step-Uke function, as it is with zero 
magnetic field. 

The electron concentration is: 



(2.10) 



^,kx 



where: 



+ 00 



d£n,,kAnfo{£)\^^,kA^\''■ (2.11) 



fo{£) is the Fermi-Dirac distribution function and 
V'ijfex (^) is the three-dimensional envelope function. 
Thus, at finite temperature , T: 



"^,fc.« = 2^-^|G,fc.(z)P 



1 



+ 00 ^ 

'^"V^ i + exp( "+^y;-^' ^' 



(2.12) 



where /x(T) is the chemical potential and is the Boltz- 
mann constant. Using Eq. (12), Eq. (10) becomes: 



(27r)2 



m*Vy = Py and m*Vz = Pz- Thus, after a little algebra 
we obtain for the accelaration and the force operators 
in the x, y and z axes, respectively: cix = +i^Vz, — 



0, Fy = and az 



1 dUjz) 



dUjz] 



. So, in the y-axis there is no force 



4-00 



da- 



va l + ea;p(^i±%V^)' 



(2.13) 



Therefore, the sheet electron concentration, is: 



(27r)2 



+ 00 

da — z 



ai+e.p r^r;T'^^ y 



(2.14) 



For a Hamiltonian like that in Eq. (1), 171*1! ~ p — qA 
22|, which in our case becomes m*Vx — Px + eBz, 



-Tin, LUUz, Ujy - 

Fz = —m*uJVx 

on the electrons, in the x-axis there is only the Lorenz 
force, while in the z-axis there is apart from the Lorenz 
force, the force due to the quantum well confining poten- 
tial. 

When there is no quantum well {U{z) = 0) the quanti- 
ties zq = "fg'j which corresponds to the z-coordinate of 
the center of the classical cyclic orbit, and xq = ^ + 
which corresponds to the x-coordinate of the center of the 
classical cyclic orbit, are constants of the motion. When 

U{z) = 0, the quantity = ^^''~^^f2,l2~^^' , which corre- 
sponds to the square of the radius of the classical cyclic 
orbit is also a constant of the motion. Thus, when U{z) 
can be ignored, electrons describe the well-known spiral 
motion. 

The algorithm used to solve self-consistently the equa- 
tions above, is divided in the following steps, {a') We 
input an initial guess for the non-magnetic potential en- 
ergy, Uin{z). (/?') We solve the envelope function Eq. 
(5) for each i and for each kx to obtain C,i^k^ {z) and 
Ei.k^ ■ Care should be taken in this step, to include all 
possible i and all possible kx which contribute to the 
electron concentration. Thus, we start with many sub- 
bands and with a wide range of k^.. This means that 
Eq. (5) must be solved many times. (7') /i(T') can be 
calculated from charge neutrality [^,^, using Eq. (14). 
{5') Thus, we can calculate, from Eq. (13) n{z) and 
ni(z) and therefore Uxc{z) [Q- (e') Now the charge 
density is known and it is used to solve the Poisson 
equation numerically [Q, to obtain Uc{z). We sup- 
pose that ^^{hulk) — 0, because there is no net charge 
in the bulk material. We take into account the differ- 
ent dielectric constants of GaAs and Al^jGai-^As [ po[ . 
Finally, we choose Uc{left bulk) — —Uq, where Uq is 
the value of the conduction band minima discontinuity. 
Uxc{bulk) — 0, because the envelope functions decay 
into the Ala;Gai_a;As barriers. Thus, U{left bulk) — 0. 
All the structures are long enough in the z-axis, so that 
bulk conditions prevail before the end of the Ala;Gai_a;As 
barriers. (^^') The output non-magnetic potential energy, 
Uout{z), can now be calculated from Eq. (2). {(') If 
Uout{z) is "very close" to Uin{z), we have finished. Oth- 
erwise, we mix Uout{z) and Uin{z) to construct the new 
Uin{z) and we return to step (/?') po|| . 

Finally, we notice that since the envelope functions 
depend on both i and kx, a quantitative calculation of 
the conductivity will involve tedious algebra, because 
the scattering matrix elements will depend on kx, too. 
This complication emerges also in the calculation of the 
screening. For this reason, although some attemps have 
already been made ||2^ , a well established transport the- 
ory for a Q2DEG with an in-plane magnetic field has not 
been developed yet. 
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III. RESULTS AND DISCUSSION 

We apply our treatment to the case of a symmetrical 
(5-doped Ala;Gai_a;As/GaAs/Al2:Gai_2;As double hetero- 
junction, in the presence of a parallel magnetic field, from 
OT up to 20T. We choose a symmetrical structure be- 
cause in this case we can observe most clearly the varia- 
tion of the electronic states induced by the magnetic field 
in the "energy dispersion picture" , in the "electron con- 
centration picture" and in the "density of states picture" . 
The structure consists of a 280A undoped Alo.25Gao.75As 
layer, a Si (5-doped Alo.25Gao.75As layer (0.28 x 10^^ 
cm^^), an undoped 250A Alo.25Gao.75As spacer, a 6OOA 
undoped GaAs well, an undoped 250A Alo.25Gao.75As 
spacer, a Si J-doped Alo.25Gao.75As layer (0.28 x 10^^ 
cm~^) and a 280A undoped Alo.25Gao.75As layer. All 
layers are assumed to have a slight unintentional accep- 
tor doping of 5 X 10^^ cm^'^. We suppose that the sam- 
ple has been illuminated and therefore all the donors are 
ionized. This is done because we want to study the ef- 
fect of the magnetic field under the condition of constant 
sheet electron concentration. Although our treatment is 
applicable to any temperature, we will apply it to T = 
4.2K. This is done because the experiments are usually 
performed at or below T = 4.2K. These material and 
structural parameters result in a sheet electron concen- 
tration, Ns= 0.54 X 10^2 cm"2_ 

First we will describe the evolution of the changes 
induced by the magnetic field to the energy dispersion 
curves, Ei{kx)- The situation is described in the lower 
parts of Fig. 1. In this particular structure, for 5 = 0, 
due to the large well width, the ground state subband 
and the first and second excited subbands are popu- 
lated, with sheet electron concentrations A^o= 0.238 x 
IO12 cm-2, Ni= 0.233 x lO^^ cm-^ and iV2=0.069 x 
10^^ cm~^, respectively. Initially, as the magnetic field 
is increased, depopulation of the higher subbands is pre- 
dicted. The second excited subband is depopulated at 
_B ~ 5T and the first excited subband at B ~ 7T. Increas- 
ing the magnetic field up to 7T, the shapes of the Ei{kx) 
dispersion curves also change. While the upper subbands 
remain almost parabolic, the first excited subband, and 
most obviously the ground state subband undergo im- 
portant changes, developing gradually local maxima at 
kx = instead of local minima aX B ~ OT. As can be 
seen from the lower parts of Fig. 1 this also happens to 
the other excited subbands for larger values of the mag- 
netic field. For B > 7T, only the ground state subband 
is populated. At this point the EQ{kx) dispersion curve 
is continuously below the chemical potential in the range 
kx = [-4 X 10®, -1-4 X 10®] m"^ This means that the 
system is still a single layer one. 

At higher magnetic fields, a transition from a single 
to a bilayer electron system occurs. This transition has 
approximately been achieved at B = 12T as can be seen 
from the lower part of Fig. Ic, but the complete separa- 
tion of the two layers is achieved at i? = 20T (see Fig. 2 



where the electron concentrations are presented). During 
this procedure, the energy separation of the unoccupied 
states (those with small \kx\) becomes hcu. This is due to 
the fact that the well width is very large and therefore 
in the central region of the well, as the magnetic confine- 
ment overcomes the well confinement, the bulk Landau 
levels dominate. This has also been predicted for square, 
analytically solvable, quantum wells when the well width 
is large enough |l6| . 

In the case of a square quantum well, the behaviour 
of the Ei{kx) curves is determined by the competition of 
the well width and the magnetic length, Ib — \/Ti/{eB) 
|l6| . When the well width is smaller than the magnetic 
length, spatial quantization dominates. The energy levels 
can be roughly classified into two types, namely, confined 
states and extended states. In this specific case the con- 
fined states in the quantum well increase parabolically 
as a function of kx |0,|l^ , while the extended states have 
an oscillating form with an "average" separation of Tiuj 
0. However, as the well width or the magnetic field is 
increased, this behavior changes. Finally, when the well 
width is larger than the magnetic length, the electron 
orbits are governed by the Lorentz force and electrons 
basically describe spiral motion. At this point the en- 
ergy dispersion curves are flat with a separation of Tiuj. 

In reference [0, the author, studing thin single quan- 
tum wells and taking as the growth axis the z-axis and 
the magnetic field in the x-axis, bypasses the dependence 
of the electronic states on the in-plane wavevector in the 
y-axis (perpendicular to B), using only ky = 0. This is 
done in order to override the large numerical cost of the 
general case. It is evident from the lower parts of Fig.l 
that such an approximation cannot be applied in our case 
because of the strong depedence of the electronic states 
on this wavevector. Moreover, for high enough values of 
the magnetic field the states with this wavevector are not 
occupied. 

The "density of states picture" is given in the upper 
parts of Fig. 1. We observe that although for B — 
IT (Fig. la) the DOS is almost step-like, there is al- 
ready a peak due to the fact that Eo{kx) has already 
developed a local maximum at kx=0, instead of the lo- 
cal minimum at i?=OT. This corresponds to a van Hove 
singularity, since '^^^l''"'^ > as we approach the critical 

point from below and '^^^^.''"^ < as we approach the 
critical point from above. The DOS of the first excited 
subband is not a "perfect step", because Ei{kx) is not 
exactly parabolic. The DOS for the second and the third 
excited subbands are "perfect steps" because E2{kx) and 
E^{kx) are parabolic. We can also see that we have three 
populated subbands. 

In the upper part of Fig. lb we present the DOS for 
B = 7T. Clearly, we can observe the depopulation of 
the first excited subband. Therefore, at this point, as 
we increase the magnetic field, the conductivity of the 
structure increases abruptly, due to the abrupt decrease 
of the DOS at the chemical potential. We also observe 
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that the total DOS is not stcp-hkc and that the second 
and the third excited subbands are not exactly parabolic. 
Moreover, since Eo{kx) and Ei{kx) have developed local 
maxima at kx=0, there are two peaks in the DOS, cor- 
responding to the two van Hove singularities. 

In the upper part of Fig. Ic we present the DOS for 
B =12T. There is a van Hove singularity at the chemical 
potential, due to the local maximum of EQ{kx) at fc^^^O. 
Therefore, at this point, as we increase the magnetic field, 
the conductivity of the structure decreases abruptly due 
to the abrupt increase of the DOS at the chemical poten- 
tial. The total DOS indicates that at the center of the 
well the bulk Landau levels start to develop. All Ei{kx) 
have already developed local maxima at kx=Q and the 
energy separation of successive subbands for small \kx\ is 
close to Tiuj. 

In the upper part of Fig. Id we present the DOS for 

B = 20T. The form of the total DOS stems from the 
combination of two factors, i.e. as we move in the en- 
ergy axis to higher energies: (a) From the two local min- 
ima of Ei{kx) up to the local maximum of Ei{kx) the 
bilayer electron system dominates, (b) From the local 
maximum of Ei{kx) up to the local minima of Ei^i{kx) 
the bulk Landau levels dominate. In this region the 
DOS has the form constant x J^ii^ - VO + ^))~^ 

= constant' x J2i,ky ^i^ - M) + f) - which is 

the DOS of a free particle in the y-axis together with 
an harmonic oscilator in the xz-plane. The energy sep- 
aration of successive subbands for small Ifc^l is equal to 

huj. 

Fig. 2 presents the variation of the electron concentra- 
tion, n{z) and of the population of the subbands, ni{z), 
as we increase the magnetic field from OT to 20T. The de- 
population of the second excited subband at i? ~ 5T and 
of the first excited subband at B ~ 7T can also be seen in 
this "electron concentration picture" . Inspection of Fig. 
2 reveals that apart from these depopulations, the form 
of n{z) changes even from IT to 7T, with a little bigger 
separation of the two parts of n{z). This separation in- 
creases with the increase of the magnetic field. At 12T 
there are still electrons in the middle of the well. The di- 
vision into two parts is complete at 20T. This means that 
the "electron concentration picture" gives a more precise 
depiction of the transision to a bilayer system than the 
"energy dispersion picture" . 

We finally give an example of an asymmetric het- 
erostructure. The structure consists of a 700A undoped 
Alo.25Gao.75As layer, a 50A Si-doped Alo.25Gao.75As 
layer (2 x 10^* cm~^), an undoped 50A Alo.25Gao.75As 
spacer, a 6OOA undoped GaAs well, an undoped 200A 
Alo.25Gao.75As spacer, a 50A Si-doped Alo.25Gao.75As 
layer (1 x 10^* cm~^) and a 6OOA undoped 
Alo.25Gao.75As layer. We suppose, again, that all the 
layers have a slight unintentional acceptor doping of 4 x 
10"'^'* cm"'' and that the sample has been illuminated so 
that all the donors are ionized. T = 4.2K. These mate- 
rial and structural parameters result in a sheet electron 



concentration, 'Ns= 1.491 x 10^^ cm'^. For B = OT, 
there are four populated subbands with sheet electron 
concentrations No= 0.742 x 10^^ cm"^, Ni= 0.406 x 
10^2 cm-2. Ar2=o.229 x lO^^ cm'^ andiV3=0.114 x lO^^ 
cm^^, respectively. 

In Fig. 3 we present the energy dispersion curves, 
Ei{kx) (lower part) and the density of states (upper 
part), for B = 5T. We notice that for this asymmet- 
ric hetcrostructure Ei{kx) 7^ Ei{—kx). The popula- 
tions of the subbands are now A^o= 1-042 x 10^^ cm~^, 
Nx= 0.324 X 10^2 cm-2, N2=Q.llQ x lO^^ cm-^ and 
A^3=0.015 x 10"'^^ cm^^, respectively. We can observe 
the "transposition and the anticrossings of the parabo- 
las" which result in a complicated form for the DOS. We 
can also see that there arc two different van Hove singu- 
larities which give the peaks in the DOS. 

Generally, both in the symmetrical and in the asym- 
metrical case, the van Hove singularities are not simply 
saddle points because the Ei{kx), as we approach the 
critical points, are not of the form —akx, a >0. The 
exact form of the dispersion curves is obtained from the 
self-consistent calculation. Anyway, as we increase the 
magnetic field, whenever the chemical potential is identi- 
fied with a van Hove singularity, the conductivity of the 
structure will decrease abruptly. On the contrary, when- 
ever there is a depopulation of a local energy dispersion 
minimum, due to the decrease of the DOS at the chemical 
potential, the conductivity will increase abruptly. 

Similar results are obtained for other values of the well 
width. It is the competition between the magnitude of 
the magnetic field and the spatial quantization - together 
with the infiuence of the number of electrons - that de- 
termines the overall behaviour of the system. Extensive 
comparison with experiment will be presented in a forth- 
coming paper. 



IV. SUMMARY 

We have calculated self-consistently the energy disper- 
sion curves, the density of states, the electron concentra- 
tion and the distribution of the electrons in the subbands, 
for Ala,Gai_a;As/GaAs/ A^Gai-^As double heterojunc- 
tions subjected to an in-plane magnetic field. 

We have systematically studied the important changes 
in the density of states, induced by the variation of the 
in-plane magnetic field. We have pointed out that these 
changes are of crucial importance for the explanation of 
the magnetoconductivity experiments. 

In the case of symmetric hcterostructures, we have 
demonstrated in the "energy dispersion picture" , in the 
"electron concentration picture" and in the "density of 
states picture" the depopulation of the higher subbands, 
the transition from a single to a bilayer electron system 
and the domination of the bulk Landau levels in the cen- 
tre the wide quantum well, as the magnetic field is con- 
tinuously increased. We have also given an example of 
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an asymmetric hctcrostructurc. 

The author wishes to thank Prof. Fabio Beltram and 
Dr. Vincenzo Piazza for many useful discussions and for 
the motivation of this work. 
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FIG. 1. Symmetric heterostructure. The energy disper- 
sion curves, Ei{kx), i=0,l,2,3 (lower parts) and the density 
of states (upper parts) drawn with a common horizontal en- 
ergy axis for (a)B=lT, (b)B=7T, (c)B=12T, and (d)B=20T. 
The DOS is in arbitrary units. The chemical potential is iden- 
tified with the zero energy. The dashed curves represent the 
DOS of each subband, ni{£), while the bold continuous curve 
represents the total DOS, n{£). 
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FIG. 2. Symmetric hetero structure. The electron concen- 
tration, n{z) (bold continuous curve) and the population of 

the subbands, nt{z) (dotted curves) for (a)B=lT, (b)B=4T, 
(c)B=6T, (d)B=7T, (e)B=12T, and (f)B=20T. 
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FIG. 3. Asymmetric heterostructure. The energy disper- 
sion curves, Ei{kx), i=0,l,2,3 (lower part) and the density of 
states (upper part) drawn with a common horizontal energy 
axis for B=5T. The DOS is in arbitrary units. The chemi- 
cal potential is identified with the zero energy. The dashed 
curves represent the DOS of each subband ni{£), while the 
bold continuous curve represents the total DOS, n{£). 
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